The Role of Friction in Compaction and Segregation of Granular Materials 
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We investigate the role of friction in compaction and segregation of granular materials by com- 
bining Edwards' thermodynamic hypothesis with a simple mechanical model and mean-field based 
geometrical calculations. Systems of single species with large friction coefficients are found to com- 
pact less. Binary mixtures of grains differing in frictional properties are found to segregate at high 
compactivities, in contrary to granular mixtures differing in size, which segregate at low compactivi- 
ties. A phase diagram for segregation vs. friction coefficients of the two species is generated. Finally, 
the characteristics of segregation are related directly to the volume fraction without the explicit use 
of the yet unclear notion of compactivity. 

PACS numbers: 45.70.Cc, 64.75.+g 



I. INTRODUCTION: 

Friction plays a significant role in packings of grains in 
a granular material, however due to the non-equilibrium 
nature of such materials their macroscopic properties do 
not trivially result from the microscopic physics of their 
constitutes. Compaction and segregation are two macro- 
scopic phenomena which occur in granular materials un- 
dergoing series of extensive operations [jj (operations of 
a statistical nature rather than "Maxwell-demon" oper- 
ations which act on single grains). Packings of identical 
hard spheres compacted by extensive operations reach 
a state of random close packing (RCP), rather than an 
ordered crystalline packing, which has a higher volume 
fraction (defined as the ratio of the sum of grain volumes 
to the volume the system occupies). Understanding the 
process leading to RCP, its geometrical properties and 
whether higher volume fractions may be obtained is of 
great importance for physics and for engineering. The 
tendency towards segregation in granular mixtures com- 
prised of grains with various mechanical properties is in- 
teresting for the physicist and disturbing for the engineer, 
for whom a homogeneous mixture is often an industrial 
need (for a review see |2j). 

Segregation in granular materials has received much 
attention in recent years (for a review see |3|]). The phe- 
nomenon is observed for grains varying in size Q , shape 
[j| , friction coefficient |(| and density @ ■ Segregation oc- 
curs due to vibration m, tapping |[|, rotation pour- 
ing |Tcj and shearing |ll| . Experiments are performed 
for mixtures of many particles of two different species 
(for a review see ^3) an d with single intruder particles 
in systems of a single species Existing theoretical 

modelin g of segregation due to rotation and to 

pouring [nj, Hg is based on kinetic phenomena: segrega- 
tion is explained as a result of different flow properties 
of the different species. Monte Carlo simulations of dy- 
namic phenomena in vibrated systems give insight into 
segregation 01 as wen as compaction dynamics of single 
species systems [Tsj . Similar dynamic phenomena have 
been captured analytically in models based on free vol- 
ume considerations 19, 20, 21]. 



This paper deals with the role of friction in compaction 
and segregation through the analysis of the static prop- 
erties of granular materials. The results obtained may 
be used in order to verify experimentally the validity of 
such models. One such proposal for the description of 
static granular materials is the analogy to the statistical 
mechanics of thermodynamic of Edwards pj , a more de- 
tailed description of which will be given in the following 
section. The central idea behind it is that even though 
the system is static and does not move with time within 
the ensemble of mechanically stable arrangement of the 
grains, we may assume ergodicity and employ statistical 
mechanics considerations for the probability of finding 
the system in any one of its states. This model requires 
the existence of an analog of temperature, referred to 
as "compactivity" (Other effective temperatures may be 
defined for jammed granular materials p2l |23||. and at- 
tempts have been done to connect them to the compactiv- 
ity |24Ll25| ). Recent experimental evidence for reversibil- 
ity in compaction processes |26j has provided justifica- 
tion for this thermodynamic analogy and have proposed 
a connection between compactivity and experimentally 
controllable quantities |27| , however the proposal remains 
controversial. 

Edwards' hypothesis is that an analog of the free en- 
ergy Y = V — XS is minimized where V is the system 
volume, S is the entropy and X is the compactivity. Our 
expectation that this formalism may predict frictional 
segregation is based on the following argument. Consider 
two systems of identical grains, such that the friction co- 
efficient of the first system, \i\ , is greater than that of the 
second system, /i2- Then S\ > S2 because every configu- 
ration available to the second system may be identically 
realized to the first system, while the converse is not true. 
This suggests that under certain circumstances a mixture 
of grains with different friction coefficients may prefer to 
segregate in order that its entropy be maximized, com- 
peting with the preference of the entropy of mixing to be 
maximized in the homogeneously mixed state. 

The statistical hypothesis has been used to investigate 
segregation in binary mixtures of species differing in size 
by mapping them to the Ising model, resulting in segre- 
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gation below some critical compactivity P, |2^| . Recent 
geometrical calculations have enabled relating the ideas 
of the statistical proposal to actual densities of granular 
systems [2^. 

In section|n]we use Edwards' statistical hypothesis to- 
gether with a simple mechanical model for the quantita- 
tive description of friction in 2D and 3D single species 
granular materials. This is combined with simple calcu- 
lations of Voronoi cell volumes for general coordination 
numbers resulting in the dependence of volume fraction 
on friction coefficient and compactivity. Section IlIII uses 
the mean-field approximation in order to describe segre- 
gation in a binary mixture of grains differing in frictional 
properties. Unlike mixtures of grains differin g in size, 
which may be mapped to the Ising model [ll |28(. fric- 
tional differences between grains result in larger entropy 
for the rougher grains. Therefore, these systems may not 
be mapped exactly onto the Ising model and segregation 
occurs above a critical compactivity and not below it. 
We then generate a phase diagram for segregation vs. 
friction coefficients of the two species. Finally, the de- 
pendence of the results on compactivity is eliminated by 
using the volume fraction as a measure for compactivity. 
Section llVl concludes with a summary of the results. 



II. SINGLE SPECIES: 
A. Statistical Model 

We use the statistical mechanics hypothesis proposed 
by Edwards for the description of jammed granular 
systems. In this formalism each mechanically stable ar- 
rangement of the grains is equivalent to a micro-state in 
statistical mechanics, and the total volume the grains oc- 
cupy plays the role of the energy. The analog of temper- 
ature is assumed to exist, is denoted X and is called the 
compactivity. We measure X in units of volume so that 
the analog of the Boltzman constant is equal to unity. 
In analogy with the canonical ensemble of states in ther- 
mal systems, the probability for the occurrence of a state 
with volume V is assumed to be proportional to e~ v l x . 

We consider a system of N identical spherical (in 3D) 
or circular (in 2D) grains. For every arrangement of the 
grains, the total volume of the system may be written as 
the sum over all grains of the Voronoi cell volume around 
each grain, Vf. 



N 



(i) 



i=l 



Average volumes of Voronoi cells, calculated neglect- 
ing spatial correlations between locations of grai ns, have 
been shown to agree with exact calculations |3(| . More- 
over, correlations have been shown to have a small ef- 
fect on the dependence of total volume on compactivity 
|29|| . Therefore, we use a mean-field approximation and 



assume the volume of every Voronoi cell is evenly dis- 
tributed between a minimal and maximal volume, v m in 
and v m ax- The geometrical and mechanical considera- 
tions determining these volumes will be presented in the 
following section, and at this stage it is only assumed 
Vmin and Vmax are uniform for all grains in the system. 

Following the analogy with the canonical ensemble, we 
may calculate the pseudo partition function: 



Z = 



-Zn/x dvi ... dv 



N 



N 



■IX . ' x -sinhf — 



N 



(2) 



where we have introduced the notations, v m id = (v m in + 
v m ax)/ '2 and At; = (v ma x - v min )/2. 

The average volume per grain may easily be derived 
from the partition function as: 



< v > = 



<V> 



N 



- N X 



,d\n(Z) 



dX 

■ X — Av ■ coth 



Av 



(3) 



This expression has been derived in and it is clearly 
seen that as X — > 0, < v >— > v m in and as X — > 00, 
< v >^ v m id- This is analogous to thermal systems, 
where at low temperatures the system is most probable 
to be found in its ground state, and as the temperature 
is increased exited states are occupied with increasing 
probability until the limit of infinite temperature, where 
all states are occupied with an equal probability, and the 
system's energy is the average energy of all these states. 



B. Mechanical Model 

We would now like to introduce mechanical and geo- 
metrical considerations to estimate v m in and v m ax , which 
must be known in order to evaluate the expression in 
eq. The minimal volume is achieved for hexago- 

nal packing in 2D and face centered cubic or hexago- 
nal close packing in 3D. The corresponding Voronoi cell 
volumes are v%® n = ^12r 2 and v^ n = y/32r 3 , where 
r is the grain radius. The resulting volume fractions 
are = nr 2 /v^ in = tt/\/12 ~ 0.91 and 

47r^,3 /„,3£> _ 



i/ v i 



tt/ 



mm 

18 ~ 0.74. 



Although purely geometric considerations determine 
Vmin, the frictional forces between the grains manifest 
themselves in v max ■ The idea behind this is that friction 
at grain contacts allows for the formation of arcs and for 
agradual decrease in the number of contacts per grain 
|31| . which in turn increases the volume of the Voronoi 
cell around every grain. 

In order to estimate the effect of friction in granular 
materials we will consider "toy-systems" consisting of a 
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small number of grains (three in 2D and four in 3D), 
calculate the effect of friction there, and use the result 
in order to obtain an approximate prediction for the de- 
pendence of v m ax on the friction coefficient, [i. This, in 
turn, may be inserted together with v m i n into eq. (jSJ) 
in order to obtain an approximate expression for the de- 
pendence of the total volume of the system, and hence 
of the volume fraction, on [i. The results will depend 
on the compactivity, whose physical significance still re- 
quires elucidation. However, even without understanding 
its significance, a few predictions for experimental results 
may be drawn from the model presented here. 

The 2D "toy-system" consists of two grains lying on 
top of a horizontal plane and a third grain lying on top 
of these (see fig. ^1 . All grains and the horizontal plane 
are assumed to be hard, and frictional forces with an 
equal coefficient of friction, /i, act at all four contacts. A 
uniform gravitational force acts downwards on all grains. 
The condition for mechanical equilibrium is that forces 
and torques acting on all grains vanish. It can easily be 
seen that this is satisfied whenever: 



sm( 



1 + cos t 



(4) 



where 8 is half the angle between contacts of the top grain 
with the bottom grains (see fig. For fi > 1 eq. Q 
is always satisfied, and any state with 7r/6<#<7r/2is 
mechanically stable. For smaller values of /i only states 
with 8 < Qmax are mechanically stable, where 8 max is 
determined from 



sin Or, 



1 + cos 6 n 



(5) 



If /i < (2 + \/3)~ 1 — 0.3, the frictional forces cannot 
hold the top grain on top of the two bottom ones, and 
these simple considerations may not be used in order to 
determine v max . In this case the volume of the Voronoi 
cell around every grain is set in the model to v m i n . Sub- 
stituting this into eq. © yields < v >= u m , n , which 
corresponds to the maximal volume fraction quoted ear- 
lier, and hence, to crystallization. Therefore, this model 
may not be used for frictionless systems, since it predicts 
crystallization at every compactivity. 

In 3D we consider three grains in an equilateral trian- 
gle lying on top of a horizontal plane and a fourth grain 
lying on top of them (see fig. As in the 2D case, all 
grains and the horizontal plane are assumed to be hard, 
frictional forces with an equal coefficient of friction, fi, 
act at all six contacts, and a uniform gravitational force 
acts downwards on all grains. Again, the condition for 
mechanical equilibrium can easily be seen to be given by 
eq. Q, however now 8 is the angle between the verti- 
cal direction and the line connecting the top grain with 
any one of the bottom grains (see fig. As in 2D, 

for small values of /i, and specifically for frictionless sys- 
tems, these considerations cannot be used to estimate the 
volume fraction, since this model predicts crystallization 




FIG. 1: The 2D "toy-system". The shaded rhombus is the 
segment of the Voronoi cell lying between two adjacent con- 
tacts. 



(a) 



(b) 



A 




FIG. 2: a) The 3D "toy-system" 
and one of the bottom grains. 



b) Section of the top grain 



regardless of compactivity, while actual 3D granular sys- 
tems do not fully crystallize, but fall into an RCP state. 

We would now like to use the results of these sim- 
ple considerations to evaluate the maximal volume of the 
Voronoi cell around every grain in a granular material. 
In such systems frictional forces allow the existence of 
large angles between contacts of every grain with its sur- 
rounding grains. We assume that the maximal angle be- 
tween contacts in a granular material depends on the 
friction coefficient in a similar manner to its dependence 
in the toy-systems. Moreover, in order to calculate the 
Voronoi cell volume around every grain we assume the 
angles between adjacent contacts of this grain with its 
surrounding grains are uniform and are given by 8 as in 
the toy-systems. Note that since we allow 8 to vary con- 
tinuously, it may not necessarily be physically possible 
to build these packings, even locally. In 2D it is pos- 
sible only if 8 is an integer fraction of 7r, while in 3D 
it is possible only if the coordination number is 4 or 6, 
which correspond to 8 ~ 1.23 and 8 ~ 0.96, respectively. 
The general expression for the Voronoi cell volume is (see 
appendix 1X1) : 



v{8) = 



nr 2 tan 6 
6 

nr 3 y/3 sin 8 tan 8 



(1+3 cos 2 fJtan-if-^/tan^tan 3 f ) 



(2D) 
(3D) 



(0) 



where a — cos 1 ((3 cos 28 + l)/4) is the angle between 
two adjacent contacts (see fig. EJ). 
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C. Results 

We combine eq. iJSJ and iJSJ to get the dependence of 
Vmax on \x. The resulting volume fraction, <I>, is now cal- 
culated using eq. © and plotted in fig. as a function 
of friction coefficient and compactivity. Compaction de- 
pends on friction only in the region 0.3 < ji < 1, where 
the mechanical model used here is relevant. As X — > 0, 
$ approaches its maximal value (determined from geo- 
metrical considerations), and as X — > oo, it approaches 
a value larger than its minimal value (determined from 
mechanical considerations) , since all possible volumes be- 
tween the minimal and maximal are equally probable. 
Obviously, <E> decreases as either /iorI increases. The 
dependence of volume fraction on friction predicted here 
may be investigated experimentally, by comparing gran- 
ular packings differing only in friction coefficient, which 
have otherwise been prepared identically, and hence may 
be reasonably assumed to have equal compactivities. 




.2 0.3 0.4 0.5 0. 

FIG. 3: Volume fraction, $, vs. friction coefficient, /i, for 
several values of the compactivity, X, for 2D (solid lines) and 
for 3D (dotted lines). Volume fractions are normalized ac- 
cording to their maximal values, $ mal = ■k/\/V2 ~ 0.91 in 
2D and $ max = ~ 0.74 in 3D. The bounding mini- 

mal and maximal volume fractions, $ m i„ and $max, are also 
plotted. 



III. TWO SPECIES: SEGREGATION 

A. Mean-Field Model 

We would now like to describe a system consisting of 
two species of grains differing only in frictional proper- 
ties. The central question we wish to address in such 
systems is whether the two species mix homogeneously, 
or segregate into separate domains. For a system of two 
species denoted A and B, we shall assume we know the 
compactivity, X, and the friction coefficients, haa, Hab 
and hbBi between two grains of type A, between a grain 
of type A and a grain of type B, and between two grains 



of type B, respectively. The considerations presented in 
the previous section may be used in order to calculate the 
maximal volumes, vaa = v max (^AA), vab = v max ((XAB) 
and vbb = v max (nBB)i of the Voronoi cells around an 
A grain surrounded by A grains, an A grain surrounded 
by B grains (or a B grain surrounded by A grains) and 
a B grain surrounded by B grains, respectively. We de- 
scribe a monodisperse system, hence the corresponding 
minimal volumes are identical for all types of grains and 
are denoted here by v m m- 

We would now like to write the partition function, Z, 
for two species and to derive from it the analog of free 
energy, Y = — X\og(Z), as a function of the concen- 
tration / = ^f, where N is the total number of grains 
and Na is the number of A grains. The constraint on 
the total number of grains of each species in the system 
causes us to view / as a local concentration which varies 
throughout the system. As in the mean-field description 
of a binary alloy, which is equivalent to the Ising model, 
a single minimum of Y(f) means the two species tend 
to get mixed homogeneously at a concentration equal to 
the global concentration determined from the number of 
grains of each species. The existence of two minima of 
Y(f) means that the system tends to separate into do- 
mains with two different concentrations. The percentage 
of the system with each of these two minimizing local con- 
centrations is determined from the global concentration 
according to the Maxwell construction (see e.g. [32]). 

The number of A grains is Na — f N and the number 
of B grains is Nb = (1 — f)N. In the mean-field approx- 
imation the number of A-A contacts is Naa — f 2 Nz/2, 
the number of A-B contacts is Nab = /(I — f)Nz and 
the number of B-B contacts is Nbb = (1 — f) 2 Nz/2, 
where z is the average number of neighbors per grain, or 
the average coordination number, which is assumed to be 
uniform for all types of grains. The contribution of every 
contact to the total volume is limited according to its 
type between 2v m i n /z and 2vij/z, with i and j denoting 
A or B for the types of the two grains in contact. The 
pseudo partition function is: 



Z 



N\ 



(fN)\((l-f)N)\ 

) 



VAA \ fN 

e~ v/x dv 



i 1 ■> j 



vab \ 2/(l-/)iV , rVBB \ (1-fVN 

e v / x dv) ( I e^dvj (7) 



Not only does Z not depend on z when z is assumed to 
be uniform for all types of grains, it can easily be seen 
that when different values of z are assigned to the dif- 
ferent types of grains, the same expression is obtained. 
This is an important result in the analogy between the 
configurational statistical mechanics of a granular mix- 
ture and the Ising model or a binary alloy. In the Ising 
model spins are arranged on an ordered lattice, all spins 
have the same number of nearest neighbors, z, and the 
total energy is determined by the states of the spins on 
all the lattice sites. In a granular system, on the other 
hand, the disordered spatial configuration determines the 



5 



volume of the system and every type of grain may have 
a different number of nearest neighbors, or a different 
coordination number, which is related to the volume of 
the Voronoi cell around it, or to the volume it occupies. 
Due to the different friction coefficients, z varies between 
the different types of grains, as is indicated in |3l| where 
the friction dependence of the average number of con- 
tacts is investigated. The number of contacts of every 
type and the contribution of every contact to the total 
volume depend inversely on z, and since only their prod- 
uct enters the total volume of the system, the resulting 
partition function does not depend on the values of z for 
the different types of grains [34( • We will soon see where 
such physical differences between the granular mixture 
and the Ising model do affect the resulting behaviour of 
the system. 

Using the partition function we now evaluate the free 
energy: 

Y(f,X,N) = -X\og(Z) = 

= XiV[/log(/) + (l-/)log(l-/) + 

+ 2/(1 - f)R(X) - log(X) - 

- fR AA (X)-(l- f)R BB (X)}, (8) 

where we have used the notations Rij(X) = 
log ( e -v min /x _ e-vtj/x) andR ( X ) = (R AA + R BB )/2- 
R AB , and the Stirling's formula has been used to evalu- 
ate logiV! for large N. Defining Avij 
that 



/ /(I _ e -A VAA /x\ 7{ Z e -Av BB /x\ \ 

R W = l-Y (l _ e -AW*) J • (9) 

We would now like to minimize Y for given overall com- 
positions, N A and N B , and for a given compactivity, X. 
Eq. (JSJ has been derived under the mean-field approxi- 
mation, therefore it describes the free energy of a region 
with uniform concentration, /. Formally we should now 
define the local concentration as a spatially dependent 
coarse grained function, /( r ), and the free energy as a 
functional of it, Y[f( r )] = J Y (/ '(~?))d~r , and to re- 
quire that f{~r) minimize Y under the constraint that 
V I f( i r)dr — N A /N. The spatial integral over all the 
system of the last three terms in eq. iJSJ, which are lin- 
ear in /, is independent of the function f(~r). Therefore 
these terms do not contribute to the minimization of Y, 
and we need only consider the contribution of the first 
three terms to Y. We thus obtain an expression similar 
to the mean-field free energy of an Ising model or of a bi- 
nary alloy |3j|. The "equilibrium" concentration is now 
determined from: 



dY 



XN 



log 



1-/ 



2(l-2f)R(X) 



which is equivalent to: 



2/-l=tanh( J R(X)(2/-l)). 



0, (10) 



(11) 



For R(X) < 1 the only solution to this equation is / = 

0. 5, which corresponds to mixing, while for R(X) > 1 
two different solutions exist and the systems segregates 
into regions with these two minimizing concentrations. 
Therefore, the condition for segregation is that R(X) > 

1, where R(X) is given by eq. ©■ 

Since R(X — > 0) = 0, no segregation occurs in the 
limit of low compactivities. Contrast this to the behav- 
ior of the Ising model and binary alloys, where phase 
separation exists at low temperatures, and specifically in 
the limit of zero temperature [221 • Here there is a min- 
imal critical compactivity, above which segregation oc- 
curs, rather than the maximal critical temperature, bel- 
low which phase separation occurs in the Ising model. 
The basic difference between the model presented here 
for granular materials and the Ising model is as follows: 
In the Ising model the energy of the system is deter- 
mined only from its topological state. Namely, once 
the topology of which element is the nearest neighbor of 
which other elements is specified, the total energy of the 
system is determined. In the mean-field approximation 
these topological states are described by the concentra- 
tion, hence the energy depends only on it. In the granular 
model presented here a topological state, specifying the 
types of contacts around every grain, allows for a range 
of volumes for the Voronoi cells around every grain and 
therefore for a range of values for the overall volume of 
the system. As in the Ising model, the topological state is 
specified by the concentration in the mean-field approxi- 
mation, however the total volume of the system depends 
both on the concentration and on the compactivity, X. 

The reason for segregation only at high compactivities 
may be understood to stem from the range of possible 
volumes in the following way. The Voronoi cells around 
the grains all have the same minimal volume, v m i n , de- 
termined from geometrical considerations. The proba- 
bility for finding a Voronoi cell in an exited state with 
a larger volume, v > v m i n , decays as e~ v l x . There- 
fore at the limit X — > all Voronoi cells are expected 
to be found in their ground state, namely to have a vol- 
ume v = v m in, independent of the friction of the grains, 
which only determines the maximal Voronoi cell volume, 
Vmax(fJ')- Hence, at low compactivities the differences 
between grains vanish and the two species mix homoge- 
neously. Segregation may occur only at high compactiv- 
ities, where exited states, which exist due to friction, 
have a greater "thermodynamic" weight. This occur- 
rence of segregation at high compactivities rather than at 
low compactivities may explain the frequent appearance 
of segregation in granular systems, which typically have 
significant compactivities, since X = corresponds to a 
crystalline packing. 

Since for R > 0, R(X) is a monotonic function (as 
can be seen from eq. @ which begins at zero, and since 
the condition for segregation is that R > 1, there is only 
a minimal critical compactivity for segregation (and no 
upper bound above which there is no segregation) so we 
may determine whether segregation is possible by inves- 
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tigating the high compactivity limit of R(X). This limit 
may easily be evaluated to be equal to 



Ra 



R(X 



oo 



( ^/Av A aAv B B , , 
V Av A B 



and the condition for the existence of a critical com- 
pactivity for segregation for given values of the friction 
coefficients is that Roo > 1. 

This condition may be understood in the following 
way: The free energy, Y = V — XS, includes a volume 
(V) term and an entropy (S) term. At high compactiv- 
ities the entropy dominates and minimizing the free en- 
ergy is equivalent to maximizing the entropy. The en- 
tropy includes two factors, a "combinatoric entropy" re- 
lated to the topological state of the system, namely which 
grain is in contact with which other grains, and a "ge- 
ometric" entropy related to the variety of volumes the 
Voronoi cell around every grain may have within one spe- 
cific topological state (This is not the case for the Ising 
model, where the topological state determines the state 
of the system, and the total entropy includes only the 
combinatoric entropy, or the entropy of mixing ) ■ 

The entropy's dependence on concentration at the limit 
X — > oo is given by 



S = -jV[/log(/) + (l-/)log(l-/) + 
+ 2/(1-/)^]. 



(13) 



The first two terms are the combinatoric entropy or en- 
tropy of mixing, while the last term is the geometric en- 
tropy. The geometric entropy plays a similar role to the 
role of the energy in the Ising model, as can be seen in 
the expression for the free energy in the Ising model (3j| : 



F = iV[fcT(/log(/) + (l-/)log(l-/)) 
- 2/(1-/) J*], 



(14) 



where k is the Boltzman constant, T is temperature, J is 
the interaction energy between neighboring spins and z is 
the number of nearest neighbors per site. Ising systems 
exhibit phase separation when kT < — Jz, which may 
now be seen to be completely analogous to the condition 
Rao > 1 for segregation in granular binary mixtures at 
high compactivities. 



B. Phase Diagram 

We would now like to generate a phase diagram for seg- 
regation and to check for which values of \iaa-, Hbb and 
Hab segregation is possible. From mechanical considera- 
tions we would expect ixab to lie between /iaa and \j,b b ■ 
Since v max (iJL) is monotonic, Ai>(/z) is monotonic as well, 
and Avab is bounded between Avaa and Avbb- First 
we notice that if hab = m o>x(fiAA, Hbb), Roo < and 
no segregation is expected. For given values of Haa and 
Hbb, Roo varies monotonicaly with [iab , and a necessary 
condition for segregation in intermediate values of [iab 



is that there is segregation at \iab 
this case 



min^/iAAi I^bb)- In 



R c 



log 



(Avaa, Avbb) \ 



i max 



min(Av A A,Av B B) I 



(15) 



Fig. 21 displays the regions in the jiAA ~ Hbb plane 
where Roo > 1, and in which segregation may occur. 
When Roo < 1 the two species are always expected to be 
mixed. When Roo > 1 segregation occurs at the limits 
of Hab = fnin(iiAAil^BB) and X — * oo. Segregation also 
occurs at larger values of \iab and at finite values of X. 
The existence of segregation may be verified for the gen- 
eral case by checking whether the expression for R(X) in 
eq. @ is greater than one. Here we only demonstrate 
our model results for the limits described above. As will 
be shown in the following subsection, in some cases segre- 
gation occurs for every compactivity. The phase diagram 
shows that mixtures of grains with close friction coeffi- 
cients will mix and what difference in friction coefficients 
is needed in order to achieve segregation. 
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FIG. 4: Phase diagram of the \iaa — I^bb plane under the 
assumption that (iab = min(fiAA, [ibb) for 2D (solid lines) 
and 3D (dotted lines) . M indicates mixing, Si and 5*2 indicate 
segregation. M regions are always mixed. In 2D S regions are 
segregated at high enough compactivity (see text). In 3D Si 
regions are always segregated, while S2 regions are segregated 
only at high enough compactivity (see section UlI Cfl . 

Within the framework of the mechanical model used 
here, frictional forces are relevant only for 0.3 < \i < 1. 
If HAA ^$ 0.3 or j-iBB ;$ 0.3 it is questionable whether 
this mechanical model is valid. If both fiAA > 1 and 
Hbb > 1 the friction has an identical effect no matter 
what the values of haa and hbb are, therefore the two 
species mix. If only one of ^aa or Vbb is greater than 
one, the maximal volume per grain of the corresponding 
species is unbounded, while for the second species it is 
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bounded, leading to Roo > 1, and hence to segregation. 
These regions are designated as well in fig. 0] 



C. Eliminating the Compactivity 

Given the mechanical properties of the grains, namely 
Haa, fJ-AB and hbb-, we have employed a mean field ap- 
proximation in order to solve the Edwards statistical me- 
chanics description of a granular mixture and to obtain 
the concentration (or concentrations in case of segrega- 
tion) the system will be found at as a function of com- 
pactivity. Even though there is evidence that it may 
be controlled experimentally [2|| |23|, the compactivity 
is still not a measurable quantity, and we would like to 
reach a description of segregation independent of the no- 
tion of compactivity. 

For a single species we have seen that the total volume 
of the system, and correspondingly its volume fraction, 
depend on compactivity. We shall now derive the com- 
pactivity dependence of the volume fraction in a binary 
mixture, and use it in order to eliminate the compactiv- 
ity from the description of segregation. This will result 
in a relation between concentrations and volume fraction, 
which may be measured experimentally. 

The average volume per grain may easily be derived 
from the partition function (eq. Q) as: 



< v > 



<V> = 1 x 2 dln(Z) 

N N dX 
X + f 2 Q AA (X) + 2f(l- 
(l-f) 2 Q BB (X), 



f)Q AB (X) 



compactivity for segregation, X Cl above which segrega- 
tion occurs and the critical volume fraction for segrega- 
tion, <1> C = Q(X C ), below which segregation occurs. An 
important point to note in the case demonstrated in fig. 
0is that in 3D $ c ~ 0.7, which is larger than the volume 
fraction of 3D RCP, <&rcp — 0.64. Granular materials 
typically do not reach volume fractions higher than that 
of RCP, therefore in this case the compactivity is always 
high enough so that the volume fraction is smaller than 
$ c and segregation always occurs. We now use this argu- 
ment in order to divide the \i AA — A 1 bb plane into regions 
where $ c > Qrcp (denoted Si in fig. which are segre- 
gated for all X, and regions where <E> C < &rcp (denoted 
S2 in fig- EJ), which are segregated only above some X c 
(or only under some <I> C ). Such considerations may not 
be used in 2D, where full crystalization is achieved ex- 
perimentally. 
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where 



Qij{X) = X 



■dRijiX) 



dX 



-v min /X 



-Vij/X 



(17) 



FIG. 5: Concentration, /, vs. volume fraction, 4>, for 
[iaa = 0.35, fj>BB = 0.75, fiAB = 0.35 for 2D (solid lines) 
and 3D (dotted lines). For low compactivities (or high vol- 
ume fractions) a single concentration is possible and the two 
species mix. For high compactivities (or low volume fractions) 
two concentrations are possible and segregation occurs. 



As for a single species, when X — » 0, Qij(X) — > u mi „ 
and < v >^ v m i n . Moreover, when X — > 00, Qij(X) 

- X, hence < v >^ v min + f 2 v AA + 2/(1 - 

f)v AB + (l- f) 2 V BB . 

The concentration / may be expressed as a function 
of X by using eq. (|11|) . Eliminating X between eq. <|11[) 
and l|16|) , we obtain / as a function of the volume frac- 
tion, $. This is plotted in fig. for typical values of 
the friction coefficients {[i AA — 0.35, [Ibb = 0.75 and 
fJ'AB = 0.35). At low compactivities the volume fraction 
is high and the system is mixed. As the compactivity is 
raised the volume fraction is reduced and at some critical 
point segregation commences, and two different concen- 
trations appear. These two values are the concentrations 
in different domains of the system. As the compactivity 
is raised above the critical point the two concentrations 
separate one from each other, enhancing the difference 
between the different domains. We identify the critical 



IV. CONCLUSIONS: 

We have used Edwards' statistical mechanics hypoth- 
esis together with a simple mechanical model to describe 
the role of friction in 2D and 3D granular materials. For 
a single species this describes the decrease in volume frac- 
tion with increasing friction coefficient and with increas- 
ing compactivity. An experimental test of the ideas pre- 
sented is most easily interpreted for systems differing in 
friction coefficient but with the same compactivity. It 
is intriguing to consider whether identical preparation 
would lead to equal compactivities for systems differing 
only in friction coefficients. 

In addition, the model has been used in order to inves- 
tigate segregation in binary mixtures of grains differing in 
frictional properties. Unlike mixtures of grains differing 



8 



in size, which may be mapped to the Ising model, fric- 
tional differences between grains result in larger entropy 
for the rougher grains; therefore, these systems may not 
be mapped exactly onto the Ising model, we find that 
segregation occurs above a critical compactivity and not 
below it. A phase diagram for segregation vs. friction co- 
efficients of the two species, which may be tested exper- 
imentally, has been generated. By eliminating the com- 
pactivity, we have also provided a relation between the 
volume fraction and the nature of mixing or segregation. 
This relation both provides an option for experimental 
validation of the model and of the statistical hypothesis 
and allows to identify mixtures which are expected to 
segregate at every compactivity. 

The geometrical and mechanical models used to de- 
scribe friction in this paper are more qualitative than 
quantitative, and more sophisticated models may be sug- 
gested. However, the qualitative results obtained here do 
not depend on their details but only on basic properties 
which all such models should have: the minimal volume 
does not depend on friction and the maximal volume and 
the number of possible states increase with friction. The 
experimental validation or invalidation of the aforemen- 
tioned results may shed light on the validity of the sta- 
tistical mechanics proposal of Edwards. 








FIG. 6: The segment of the 3D Voronoi cell around the grain 
centered at O and blocked between the three contacts A, B 
and C. The tangent planes at the contact points meet the 
planes formed by the directions to the contact points at D, E 
and F and all meet at G. 
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APPENDIX A: VORONOI CELL VOLUMES 

This appendix describes the calculations leading from 
the mechanical models in 2D and in 3D to the volume of 
the corresponding Voronoi cells around the grains, given 
in eq. @ . In 2D the Voronoi cell around a grain is formed 
by the tangents to it at the contact points with the sur- 
rounding grains. The segment of the Voronoi cell lying 
between two contacts is the shaded rhombus in fig. ^ 
which has an area of r 2 tan 9, where r is the grain radius 
and 9 is half the angle between two adjacent contacts. 
If all angles between adjacent contacts are equal to 29, 
the Voronoi cell is comprised oiir/9 such segments. Even 
though this is an integer number only for discrete values 
of 6, we use this for continuous values of 9. The resulting 
area of the Voronoi cell is = 7rr ^ an9 . 

In 3D we use a triangulation of the contact points on 
the surface of every grain, and in analogy to the 2D case, 
we consider the segment of the Voronoi cell bounded be- 
tween three contacts (see fig. . We assume the three 
contact points A, B and C form an equilateral triangle, 



therefore ZAOG = ZBOG = ZCOG = 9, where O is the 
center of the grain and G lies on the line directed from 
O to the center of the triangle A ABC . The segment of 
the Voronoi cell is the region bounded between the three 
directions OA, OB and OC (or the planes AOB, BOC 
and CO A) and the tangent planes ADGF, BEGD and 
CFGE (which due to symmetry meet at the point G, 
which is located along the direction from O to the center 
of the triangle A ABC). The segment of the Voronoi cell 
may be divided into six tetrahedra of equal volume all 
having the common edge OG: OADG, OAFG, OBDG, 
OBEG, OCEG and OCFG. Its volume is hence given 
by 



Jseg 



6vq 



CEG 



r 2 tan 9 ■ 



c SocgHe 

\/3r sin 9 
1 + 3 cos 2 9 ' 



(Al) 



where Soc'G is the area of the triangle A OCG and He 
is the distance of E from the plane of A OCG. In order 
to calculate the total volume of the Voronoi cell we cal- 
culate the number of segments it is comprised of through 
the solid angle every such segment occupies. The an- 
gle ZAOB = ZBOC = ZCOA = a may be calculated 
from the scalar product of any two of the vectors OA, 
OB and OC and is given by cos a = (3 cos 29 + l)/4. 
The solid angle formed by these three vectors is ft — 

4tan _1 ^ v /tan(3a/4) tan 3 (a/4)^) . Therefore the Voronoi 
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cell volume is 

3D _ 4tt 

v vor — v seg~Q 



;rr J v / 3sin^tan0 



" -(A2) 



(1 + 3 cos 2 6) tan" 1 (yW^tan 



3 a 
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